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The formation of ultralong-range Rydberg molecules is a result of the attractive interaction be¬ 
tween Rydberg electron and polarizable ground state atom in an ultracold gas. In the nondegenerate 
case the backaction of the polarizable atom on the electronic orbital is minimal. Here we demon¬ 
strate, how controlled degeneracy of the respective electronic orbitals maximizes this backaction 
and leads to stronger binding energies and lower symmetry of the bound dimers. Consequently, the 
Rydberg orbitals hybridize due to the molecular bond. 


The geometrical structure of molecules determines 
their physical and chemical properties. The shape of an 
individual molecular bond can be explained by the con¬ 
cept of hybridization. The mixing of nearly degenerate 
atomic orbitals leads to a new hybrid orbital, which al¬ 
lows to analyze structures of many basic molecules like 
carbon dioxide, ammonia or water. One of the first 
molecules described this way is methane [1]. In CH 4 the 
carbon electrons, initially in s or p configuration, are re¬ 
arranged into four hybrid orbitals sp^ while creating the 
molecular bond. The new hybrid orbital contains 25% of 
s character and 75% of p character. 

Here, we present the hybridization of a Rydberg elec¬ 
tron orbital induced by the formation of an ultralong- 
range Rydberg molecule. Rydberg molecules consist of 
a single Rydberg atom bound to one or more ground 
state atoms. The bond in this type of molecules origi¬ 
nates from the elastic scattering between the slow Ryd¬ 
berg electron and a ground state atom. Diatomic Ryd¬ 
berg molecules have been observed for S-states nig, D- 
states H [5] and P-states 0 [g in gases of Rb or Cs. 
Polyatomic Rydberg molecules have been recently pho- 
toassociated from a gas of Rb 0. The presence of the 
neutral perturber inside the Rydberg electron orbital can 
cause admixing of a nearly degenerate hydrogenic mani¬ 
fold and energetically close by /-states at the location of 
the ground state atom 0. As a result, these so-called 
trilobite molecules, can possess a giant permanent elec¬ 
tric dipole moment. For Cs, which has an almost inte¬ 
ger quantum defect for the S-state, this admixture can 
amount to 90% and lead to a permanent dipole moment 
of kilo-Debye, which was shown in m for principal quan¬ 
tum numbers n=37, 39, 40. In Rb, the hydrogenic man¬ 
ifold is energetically much further away. Therefore the 
hydrogenic admixture was less then 1 % and the observed 
permanent dipole moments were on the order of 1 Debye 
for n=35, 43 m- In this paper, unlike in the experi¬ 
ments with trilobite states, we do not mix many states 
with different principal and angular quantum numbers. 
We demonstrate the perturber induced formation of a 
new hybrid orbital, which is a linear combination of only 
two mj-states with known orbital shapes. For this rea¬ 
son, it makes it a textbook-like example of tunable hy¬ 


bridization in the field of Rydberg molecules. The ratio 
of the contributing states can be precisely controlled by 
an applied electric field. This allows for tuning the con¬ 
tribution of one mj-state to the hybrid orbital from 0 to 
100% over the crossing of the two m^-states. 

In the absence of an external field, the magnetic sub¬ 
states of the Rydberg energy levels are degenerate. This 
degeneracy can be lifted by applying a magnetic field 
B. The initially degenerate states splits then into 2j + I 
different Zeeman levels labeled by the projection of the 
total angular momentum on the quantization axis rrij. 
Individual Wj-states can be brought back to degeneracy 
by applying an additional electric field S aligned paral¬ 
lel to the B-field. In Fig. calculated Stark maps in 
the vicinity of the 42D state in an external magnetic 
field of B=13.55 G are shown. The eigenenergies are ob¬ 
tained by diagonalizing the full Hamiltonian consisting of 
an unperturbed atomic Hamiltonian and two terms ac¬ 
counting for the interaction with the magnetic and the 
electric field. The splitting between states of different 
quantum numbers /, with / < 3 is caused by their quan¬ 
tum defects nang. These states exhibit a quadratic 
Stark effect, as can be seen for the 42 D 5/2 state, which 
splits into six rrij levels, shown in the magnification of 
Fig. 0 States of the manifold with / > 3 show a linear 
Stark effect, which results in a fan-like Stark structure. 
The magnitude of the corresponding shifts depends on 
the absolute value of rrij. Since the electric and mag¬ 
netic dipole operators do not couple states of Arrij ^ 0 
for E ||B, the mj=5/2 and TOj=I/ 2 states become de¬ 
generate at i5=I.I35V/cm exhibiting real crossing. The 
absence of coupling between the atomic states provides 
a clean two level system, suitable for the investigation of 
hybridization of the electronic orbitals due to perturber 
induced coupling. 

To describe the interaction of a Rydberg atom with a 
ground state atom in an external homogeneous electric 
and magnetic field, we write down the Hamiltonian: 

p2 

Htot = Hq -|- Hb + Hg + — -I- I4(r, R), (I) 

T~t 2 

where ^ represents the kinetic energy of the ground 
state perturber in a center of mass frame, where M is 
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FIG. 1: Calculated Stark map in the vicinity of the 42D state for ®^Rb. The electric field is aligned parallel to the magnetic 
field 5=13.55 G. All energies are plotted relative to the energy of the 42D^/2,mj = 1/2 state at £=0. The states are labeled 
with quantum numbers of the states they adiabatically unite with at £=0. In the magnification the sixfold splitting of the 
42 D 5/2 is shown. The green circle marks the investigated crossing of the mj = \/2 (blue) and mj=5/2 (red) energy levels. The 
part of the spherical harmonics that is relevant for the molecule formation at £—0 is depicted close to the corresponding levels. 
The shaded areas a. (magenta) and b. (light blue) correspond to the measurement ranges of Fig. 


the reduced mass, Hq is a field-free electron Hamiltonian 
and Hb and Hs account for the interaction with the elec¬ 
tric and magnetic field, respectively. In this approach the 
total three body system is described by two relative po¬ 
sitions R of the ground state atom and r of the electron, 
taken with respect to the Rydberg ionic core as well as the 
respective momenta P and p. Additionally, particles are 
treated as point-like. The pseudopotential I4(r,R) de¬ 
scribing the interaction between the low-energy electron 
and the perturbing atom takes, in the s-wave regime, the 
form m 

K(r,R) = ?^^5(r-R), (2) 

nie 

where is the s-wave scattering length and m-e is the 
electron mass. For ®^Rb the triplet scattering length is 
negative and as a consequence the potential (|^ is at¬ 
tractive. Taking into account the local electron den¬ 
sity, the resulting potential is proportional to |'1>(R)|^ 
[S]. Its modification due to the p-wave shape reso¬ 
nance can result in a binding by internal quantum reflec¬ 
tion m and in the creation of butterfly-type molecules 
m- Momentum-dependent corrections to the scattering 
length can be calculated using a semiclassical approxi¬ 
mation m- We apply a Born-Oppenheimer approxima¬ 
tion and express the total wavefunction as 'I'(r,R) = 
V^(r,R)</.(R), where '0(r,R) is the electronic molecu¬ 
lar wavefunction in the presence of the perturber and 
^(R) describes the molecular rovibrational state. This 
allows for computing the adiabatic potential energy sur¬ 
faces (APES) for the fine structure states. Subsequently, 
one can calculate the molecular binding energies and 
the rovibrational molecular wavefunctions by solving 


the Schrodinger equation using the previously calculated 
APES. This method was used before for theoretical stud¬ 
ies of molecules in external magnetic and electric fields 
[TBIED] and was experimentally verified in [J]. 

Here, we are interested in the mixing of two orbitals 
mj=5/2 and mj=l/2 due to the presence of the neutral 
atom inside the Rydberg electron orbit, which manifests 
itself in a change of the molecular binding energy E^. We 
treat the molecular potential I4(r, R) as a perturbation 
in first order to the Hamiltonian H' = Hq + Hb + H£. 
Eor the calculation of the full Hamiltonian we use the 
basis of 42 D 5 / 2 , raj=hl2 and 42 D 5 / 2 , mj=ll2 electron 
wavefunctions 


?/l/2(r) = Ra 2 d { y ) Ei/2(0) = ^R42D{r) T2°(0) (-3) 

V’5/2(r) = i?42r>(r)T5/2(0,«<>) = i?42D(r) 1^(6, </), 

where we include the respective Clebsch-Gordan coeffi¬ 
cients and the spherical harmonics (shown in Eig. in 
accordance to triplet state of the electron involved. The 
radial wavefunctions are the same for both states. Since 
we take the unperturbed wavefunctions as a basis, there 
is no dependence on R left. Consequently, we obtain the 
full Hamiltonian 


/l^l/2p+^l/^ ^1*/2^5/2 \ 

y ^i/ 2^5*/2 1^5/2^ + A2/ryy 


(4) 


where 77 = 7?(R) = / drI4(r,R) |i?42D(r)P and Ai and 
A 2 are the summed-up energy shifts arising from H'. 
By diagonalizing the Hamiltonian Q we obtain the new 
eigenenergies of the system. The degeneracy between the 
basis states occurs for Ai = A 2 . In the experiment we 
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can control the energy difference Ai = A 2 — A 2 and thus 
the states mixing by tuning the electric field. 

We excite Rydberg atoms in a two-photon process from 
a magnetically trapped spin-polarized ultracold cloud of 
®^Rb atoms in the 55'i/2,.F = = 2 state [5T] with 

an applied magnetic offset field of B=13.55G. The typ¬ 
ical temperature of the cloud is 2 /iK and the peak den¬ 
sity is on the order of 10^^ cm“^. The 780nm laser drives 
the lower transition. The detuning from the intermedi¬ 
ate 5 P 3/2 state is 500 MHz. The 480nm laser driving 
the upper transition is on constantly. After each 50 /iS 
long 780nm laser pulse we ionize Rydberg atoms and 
molecules with an electric field and detect the ions on 
a microchannel plate detector. We perform 800 cycles 
of excitation, ionization and detection in one cloud while 
scanning the frequency of the red laser. This way we ob¬ 
tain one full spectrum of the Rydberg signal for a given 
value of the electric field. We adjust the polarization of 
the light driving the upper transition such that the inten¬ 
sities of the two atomic lines close to the crossing (Fig.[^ 
are comparable. 

In Fig. we show the measured Stark maps in the 
vicinity of the mj=5/2 and mj=l/2 atomic states’ cross¬ 
ing, marked in Fig. by the green circle. Due to stray 
electric fields present in the experiment, the electric and 
magnetic field are not perfectly parallel. The angle be¬ 
tween E and B is on the order of a few degrees, which in 
principle results in a very small coupling between the two 
atomic states. Nevertheless, given the resolution of our 
experiment, the considered atomic states appear as a true 
crossing. The molecular states are shifted to the red with 
respect to the atomic line. This energy difference corre¬ 
sponds to the molecular binding energy E^. Therefore, 
we can identify the molecular states by comparing their 
binding energies with the calculated values |4] from trac¬ 
ing them back to the zero electric field (Fig. Ep)- The 
molecular state visible in the Fig. originates from a 
ground state toroidal molecule bound in the equatorial 
plane (0 = 7r/2) of the mj=l/2 scattering potential (see 
inset of Fig. with Eb=3.6MHz at £—0. The molec¬ 
ular ground state bound in the axial lobe (0 = tt, 0) of 
the same scattering potential occurs at Eb=13.7MHz. 
This axial molecule follows in parallel the trace of the 
mj=l/2 atomic line up to the degeneracy point, whereas 
the toroidal bends around the crossing and asymptoti¬ 
cally follows the mj=5/2 atomic line. In the experiment 
we are able to address the individual rovibrational molec¬ 
ular states visible in the Fig. However, since the 
excited toroidal states are blurred around the atomic de¬ 
generacy, we focus our analysis on the molecular ground 
states. 

Assuming no hybridization, the molecular states of the 
mj=I/2 and mj=5/2 would occur at a constant bind¬ 
ing energy Eb(i£’= 0) with respect to the corresponding 
atomic lines regardless of the applied electric field. In 
Fig. [2] such unperturbed ground state toroidal molecule 


traces, plotted in reference to the calculated atomic en¬ 
ergy levels, are shown for the mj=I/2 and mj=5/2 states 
with continuous brown and magenta lines, respectively. 
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FIG. 2: Stark map in the vicinity of (a.) the degeneracy 
of the two atomic lines mj=5/2 and mj = l/2 and (b.) the 
mj=l/2 state followed from the zero electric field up to the 
crossing. Computed atomic states positions, shown before in 
Fig.g are depicted with continuous white lines. All theoret¬ 
ical molecular lines are plotted with respect to these atomic 
energy levels. The hypothetical positions of the ground state 
toroidal molecular states, with no hybridization of the elec¬ 
tron orbital, are indicated with brown and magenta lines, for 
mj=5/2 and mj = l/2, respectively. The binding energy of the 
toroidal molecule belonging to the mj = l/2 state changes in 
the vicinity of the atomic crossing. Calculated positions of 
the ground state toroidal and axial molecules are indicated, 
for the sake of visibility, with red and pink dotted lines, re¬ 
ceptively. In addition excited rovibrational toroidal states are 
visible in b. However, it is not possible to accurately trace 
them for higher electric fields. We identify the weak diago¬ 
nal lines visible at higher electric field as laser sidebands with 
no physical meaning. Note that the atomic lines are satu¬ 
rated due to the chosen color scale. The discontinuity in the 
mj=l/2 atomic line above 1.2V/cm is an experimental arti¬ 
fact. 
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These two states, described by the wavefunctions (|^, 
form a basis for our two level model. Consequently to 
predict the behavior of the perturbed hybrid states, we 
calculate the eigenenergies E± of the Hamiltonian 
For a given molecular state, the obtained analytical solu¬ 
tion depends only on the energy difference A between the 
mj=l /2 and mj=5/2 states and a parameter 77 , which is 
determined from the fitting of E+ to the experimental 
data at £=0. In the experiment the energy difference 
A around the crossing depends to a good approximation 
linearly on the electric field. E_, which would correspond 
to the toroidal molecule of the mj =5/2 state was not ob¬ 
served in this experiment. By calculating E_|_, i. e. the 
binding energy, as a function of A(£’), we reproduce the 
behavior of the toroidal and axial molecular lines (Fig. 

with high accuracy. Note that the modeled molecular 
traces are plotted with respect to the computed atomic 
positions. The scattering potential couples the unper¬ 
turbed toroidal molecular lines, although the respective 
atomic states cross without any level repulsion between 
them. This leads to orbital mixing and therefore a change 
in the Rydberg electron probability density, which we 
observe as an increased binding energy of the toroidal 
molecule. In case of the axial molecule, the off diago¬ 
nal coupling terms in the Hamiltonian vanish. This 
is due to the fact that the torus-like m_,=5/2 spherical 
harmonic is zero along the axial direction and thus does 
not modify the mj=l/2 state in this direction. For this 
reason the binding energy stays constant. 

In Fig.j^the measured and calculated (E+) binding en¬ 
ergies of the considered molecules versus the electric field 
are shown. Here the binding energy is determined as the 
difference between the measured atomic and molecular 
energy. The binding energy of the toroidal molecules in¬ 
creases from Eb(£=0 ) of the mj=l/2 state, reaches its 
maximum of Eb=11.7MHz and asymptotically decreases 
to Eb(£^= 0) of the mj=5/2 state. This transition can be 
attributed to an increase of the admixture of the mj=5/2 
state from 0 to 100%. The maximum binding energy cor¬ 
responds to an orbital consisting of 50% mj=l/2 char¬ 
acter and 50% mj=5/2 character. The resulting shape 
of the orbital is shown in the inset of Fig.[^ The hy¬ 
bridization introduces a ^-dependence to the Rydberg 
orbital, even though the asymptotic states were spheri¬ 
cally symmetric along the z-axis. For the unperturbed 
axial molecule, the shape of the Rydberg orbital remains 
the same. 

Directly at the crossing the experimental binding en¬ 
ergy of the toroidal molecule and the calculated one de¬ 
viate by 9%. The predominant reason of the discrepancy 
around the crossing, for £ > 0.8V/cm, are difficulties in 
determining the atomic peak position caused by its non- 
Gaussian shape. Additionally, for £" > 1.2 V/cm fitting of 
the molecular state position becomes problematic. More¬ 
over, for the calculation we only take the unperturbed 
states and not the real Stark states into account. Finally, 
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FIG. 3: Measured binding energies of the axial (orange points) 
and toroidal ground state molecules (light blue and violet 
points) versus the electric field. The binding energies are 
extracted from two data sets shown in Fig. (a. and b.) 
and hence indicated with two colors for the toroidal molecule 
data. A ground state atom inside the Rydberg electron in 
the toroidal plane mixes two orbitals. As a result the bind¬ 
ing energy of the mj =l/2 toroidal molecule increases, reaches 
its maximum directly at the crossing and then asymptotically 
decreases to the Eb value of mj =5/2 toroidal molecule. The 
relevant part for molecule formation of the spherical harmonic 
of the hybridized Rydberg orbital at the crossing is shown in 
the inset together with the two asymptotic shapes. The error 
bars are determined as two standard deviations of the fit. 


more levels could be included into the analysis. In view 
of these facts, the agreement between the experimental 
data and our simple model is remarkable. A more com¬ 
plicated analysis would not lead to a significantly better 
agreement due to the purely experimental constrains. 

We have shown the hybridization of the Rydberg elec¬ 
tron orbital due to the molecule formation around the 
crossing of the atomic lines mj=I/2 and m^=5/2. The 
backaction of the bound perturber increases the bind¬ 
ing energy and consequently changes the shape of the 
electron orbital. Mixing more than two electron orbitals 
of known shapes by bringing them to degeneracy could 
result in even more complex asymmetric electron config¬ 
urations. Furthermore, the influence of a different type 
of perturber, like an atom of another element or a het- 
eronuclear molecule [52], as well as the effect of a few 
perturbers on the electron orbital could be investigated. 
Rydberg electron orbital hybridized due to thousands of 
the ground state perturbers could be observed as an im¬ 
print on a Bose Einstein condensate [23]. Finally, en¬ 
gineering a state by placing the perturber precisely at 
a desired position might also be feasible in experiments 
with individual atoms in microtraps. 
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